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Abstract 

The CRESST experiment monitors 300 g CaW04 crystals as targets for particle interactions in an ultra 
low background environment. In this paper, we analyze the background spectra that are recorded by three 
detectors over many weeks of data taking. Understanding these spectra is mandatory if one wants to further 
reduce the background level, and allows us to cross-check the calibration of the detectors. We identify a 
variety of sources, such as intrinsic contaminations due to primordial radioisotopes and cosmogenic activation 
of the target material. In particular, we detect a 3.6 keV X-ray line from the decay of 41 Ca with an activity 
of (26 ± 4) ^Bq, corresponding to a ratio 41 Ca/ 40 Ca = (2.2 ± 0.3) x 10~ 16 . 
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1. Introduction 

The Cryogenic Rare Event Search with Super- 
conducting Thermometers CRESST [1, 2] aims to 
detect rare nuclear recoils from elastic scattering of 
Dark Matter particles [3]. Since less than one such 
recoil is expected per kilogram of target mass and 
week of exposure, efficient shielding of the target 
against ambient radioactivity is mandatory. 

To this end, the CRESST experiment is located 
in the Laboratori Nazionali del Gran Sasso under 
an average of 1400 m of rock overburden, where the 
cosmic muon flux is reduced by about six orders 
of magnitude with respect to the surface. Addi- 
tionally, the target is surrounded (from the outside 
to the inside) by 45 cm of polyethylene to reduce 
the neutron flux, a radon box that is constantly 
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flushed with nitrogen gas to reduce the concentra- 
tion of radioactive radon and its daughters inside 
the shielding, as well as 20 cm of lead and 14 cm 
of pure copper against the external gamma back- 
ground [2]. The understanding of the remaining 
radioactive background present in the experiment 
demonstrates its excellent performance and is a nec- 
essary step to improve the sensitivity of CRESST 
to rare processes. 

2. Setup 

The target consists of individual 300 g CaW04 
crystals that are operated as calorimeters [1]. Up 
to 33 such detectors can be arranged in a com- 
pact layout [2]. To enhance the temperature signal 
following an interaction, the crystals are cooled to 
~ 15 mK where heat capacities are low. The high 
frequency phonons that are created in interactions 

May 26, 2009 



in the crystal are collected by a thermometer, which 
is a thin tungsten film that is evaporated directly 
onto the crystals. This film becomes superconduct- 
ing at these temperatures and is stabilized in its 
transition to the superconducting state. A small 
temperature rise results in an increase of resistance 
that is measured with a SQUID-based readout cir- 
cuit [4]. 

A resistive heater structure is evaporated onto 
each thermometer. This structure is used to pro- 
vide a constant heating bias as well as for periodic 
injection of additional heat pulses. Each such in- 
jected heat pulse briefly drives the thermometer 
out of its transition. The amplitude response of 
the thermometer to the injected heat is controlled 
to stay at a constant value with small adjustments 
to the heating bias by means of a feedback system. 
In this way the operating point of the thermometer 
is stabilized over several months to within the re- 
quired precision [1, 2]. This can also be seen from 
the sharpness of gamma lines discussed here. 

3. Pulse Height Evaluation 

Our thermometers yield rather slow pulses with 
rise times of about 1 ms. These pulses are digitized 
with a time base of 40 fis. 4096 such samples com- 
prise a record, of which the first 1024 samples are 
taken from the time before the trigger to define the 
baseline of the pulse. 

Template pulses are fitted to the record in order 
to evaluate the amplitude of the pulse as a measure 
of the deposited energy. These templates are an 
average of appropriate pulses, often obtained from 
a gamma calibration line. Since the SQUIDs only 
measure relative changes in current, the baseline 
needs to be defined individually for each recorded 
pulse. Hence, the free parameters in the fit are the 
level of the baseline, the onset of the pulse, and its 
amplitude. 

The pulse samples the transition curve of the 
thermometer from its superconducting to the nor- 
mal conducting state. For interactions with recoil 
energies below ~ 200 keV, the resultant pulse gener- 
ally samples only the linear region of the thermome- 
ter's transition curve between the superconducting 
and normal conducting state. For interactions with 
larger energies, the pulse samples a region of the 
transition curve that is non-linear, distorting the 
resulting pulse shape. Therefore, in the template 
pulse fit, only those parts of the pulse are included 



that show a linear behavior. This linearizes the en- 
ergy response of the fitted amplitudes over a wide 
energy range, even up to the MeV range [5]. 

4. Energy Calibration 

To reach the cryogenic temperatures, the detec- 
tors are enclosed by a fivefold thermal shield made 
from copper. Calibration sources are placed outside 
this shield, which has a total thickness of 12 mm of 
copper. This requires gamma energies > 100 keV 
in order to penetrate this barrier. Hence, a method 
is required to take the energy calibration to lower 
energies. 

To this end, additional heater pulses correspond- 
ing to a number of fixed energies are injected into 
the heater structure every 30 seconds. The pulse 
height of the response pulse from the thermometer 
is also evaluated with a template fit. Non-linearities 
are due to the shape of the transition curve and 
due to the readout circuit. They manifest them- 
selves in the amplitude of the response pulse, but, 
for given amplitude, not in its shape. Therefore, 
the response pulses give the correct energy assess- 
ment even if there are small differences in the pulse 
shape of heater response pulse and particle pulse. 

The heater pulses are calibrated at relatively high 
energies, typically with the 122 keV or 136 keV line 
from a 57 Co source. The amplitude response of the 
thermometer as a function of energy injected into 
the heater can then be fitted with a low-order poly- 
nomial function. Figure 1 depicts this procedure 
for the calibration of particle pulses in the com- 
plete energy range covered by the heater pulses. In 
the following we identify a variety of lines at the 
expected energies, thus validating our calibration 
procedure. 

5. Calibration Spectrum 

Figure 2 shows a typical spectrum from a 57 Co 
calibration. The two strong peaks at 122 keV and 
136 keV originate from gamma transitions in the 
57 Co source. They show pronounced shoulders to- 
wards lower energies, which originate from photons 
that are Compton scattered in the copper thermal 
shields and the detector support structures that 
surround the detectors. 

A gamma can interact in the crystal via the pho- 
toelectric effect, ejecting an electron from its shell. 
The X-ray produced in the subsequent recombina- 
tion may escape from the crystal. In this case, we 



10 4 



|llll|llll|llll|llll|llll|llll|llll|llll|im|llll|llll|llll|llll|llll|llll|lll± 




Figure 2: Energy spectrum from a high statistics 57 Co calibration (detector Daisy, run 31). 
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Figure 1: Calibration at low energies with heater pulses of 
various discrete energies. Fitting a low order polynomial 
to the thermometer response (small points, including error 
bars) gives the transfer function (solid line) between injected 
voltage (scale on the left) and fitted amplitude (scale on 
the bottom). The 136 keV peak gives a calibration point 
(dotted line) for the injected energy (scale on the right). 
Then, for each particle pulse, the particle's energy is found 
by evaluating the transfer function at the fitted amplitude 
(dashed line). 



measure the energy of the original gamma minus 
the energy of the X-ray (escape peak) . If the escap- 
ing X-ray reaches a neighboring crystal, its energy 
will be measured there. In figure 2 we observe a set 
of peaks in the energy range between 50 keV and 
80keV, originating from the dominant K a x, K a2 , 
Kpi/Kps and escape lines of tungsten X-rays. 
The origin of the various observed lines are labeled 
in the figure. In addition, a line is visible at 8keV 
from copper fluorescence due to interactions in the 
copper surrounding the cryostat (see section 6.6). 



The calibration spectrum allows confirmation of 
our energy calibration method. This is shown in 
table 1, where measured energies are compared to 
the literature values [6] for a few lines. The values 
agree within the resolution of our detectors. 



line 




literature 


reconstructed 


resolution 


origin 




value 


energy 


(la) 


Cu K a 




8.0 keV 


8.3 keV 


0.1 keV 


122 keV - 


K/31 


54.9 keV 


55.0 keV 


0.4 keV 


WK Ql 


59.3 keV 


59.6 keV 


0.5 keV 


122 keV- 


K Ql 


62.8 keV 


63.0 keV 


0.3 keV 


57 Co 




122.1 keV 


122.2 keV 


0.6 keV 


57 Co 




136.5 keV 


136.5 keV 


0.6 keV 



Table 1: Literature values, reconstructed energies and reso- 
lution of a few lines observed in the cobalt calibration. All 
values agree and validate our energy calibration procedure. 
The energy resolution is seen to be better than 1 keV. 



6. Background Spectra 

In the following, spectra from exposure to back- 
ground radiation alone are analyzed, from one de- 
tector during the prototyping phase of the experi- 
ment [1] (figures 3, 4) and two detectors during the 
commissioning phase [2] (figures 5, 8 and 6). All 
these crystals were grown by the General Physics 
Institute, Moscow, Russia, but in different batches. 

6.1. Strontium and Yttrium 

Strontium belongs to the same chemical group as 
calcium and can thus easily be incorporated in our 



CaWC>4 crystals in the Czochralski growth process. 
Indeed this is the case, as can be seen in figure 3 
where we show the spectrum of one crystal up to 
1 MeV. We observe a continuous background which 
can be attributed to the beta decay of 90 Sr (with 
an endpoint energy of 546 keV) together with the 
subsequent beta decay of 90 Y (2282 keV endpoint 
energy). 
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Figure 3: Background spectrum observed with one crystal 
(m = 306.8 g) from an exposure to background radiation 
of 15.00 kg d early 2004 (detector Daisy, run 28). The black 
line is a fit above 80 keV of a paramctrization of the 90 Sr/ 90 Y 
beta spectrum to the data. 

A fit above 80 keV to a parametrization of the ex- 
pected beta spectrum from 90 Sr in equilibrium with 
90 Y [7] can be seen to describe the data in this en- 
ergy range. From the fit we deduce an activity of 
(4.37 ± 0.15) mBq of 90 Sr for this crystal. In ad- 
dition to this common and continuous background, 
the three crystals show quite different features in 
the energy range below 100 keV. 
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Figure 4: Zoom on the energy region below 100 keV of the 
spectrum shown in figure 3. 

6.2. Lead 

Isotopes from natural decay chains are present 
in the crystals and can be identified by their alpha 
decays [5]. A step in the natural decay of 238 U is 
210 Pb, which beta-decays into 210 Bi with an end- 
point energy of 63.5 keV. In 84% of these decays 
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Figure 5: Background spectrum of a second crystal (m = 
305.5 g) from an exposure of 21.20 kg d in 2007 (detector 
Verena, run 30). 
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Figure 6: Background spectrum recorded in 2007 from a 
third detector (m = 307.0 g) with an exposure of 22.81 kg d 
(detector ZORA, run 30). 

a 46.5 keV gamma is emitted. If the decays hap- 
pen in the vicinity of our crystals, we observe a line 
at 46.5 keV (see figure 4), which reflects an exter- 
nal activity of (31 ± 5) /iBq. If, on the other hand, 
210 Pb is a contamination intrinsic to the crystal, our 
calorimetric measurement gives the energy of the 
gamma plus that of the emitted electron. Hence, 
we observe a beta spectrum starting at 46.5 keV 
that extends up to 63.5 keV, as can prominently be 
seen in figure 5. The activity of this decay due to 
sources both internal and external to the crystal is 
much higher, namely (672 ± 28) /iBq. Figure 6 also 
shows this lead feature for the third crystal with an 
activity of (241 ± 25) //Bq. 

6.3. Actinium 

A similar case is that of 227 Ac, a step in the nat- 
ural decay of 235 U. 227 Ac beta-decays with an end- 
point energy of 44.8 keV into Th, where two ex- 
cited levels may lead to the emission of 24.5 keV and 
9.3 keV gammas. This results in two beta spectra, 
one starting at 9.3 keV, and one at 24.5 keV, as can 
prominently be seen in figure 6 but also in figure 5. 
This contamination is not observable in the spec- 
trum shown in figure 4, which indicates that this 



detector contains less impurities from the decay of 

235jJ 

6.4- Activated Tungsten 

If the raw material or the crystals themselves are 
exposed to cosmic radiation, the tungsten in the 
CaWC>4 crystals can be activated. A possible chan- 
nel is 182 W(p, a) 179 Ta. 179 Ta then decays via elec- 
tron capture with a half-life of 1.8 years into 179 Hf. 
The energy signature of this decay in our crystals is 
the binding energy of the captured electron, mostly 
from the K-shell, which is i?K.Hf = 65.4 keV. The 
activation via 183 W(p, t) results in 181 W. This de- 
cays with a half-life of 121 days via electron capture 
into an excited 181 Ta nucleus, which de-excites with 
emission of a 6.2 keV gamma. Hence, the dominant 
energy signature of this decay is the binding energy 
of a tantalum K-shell electron, Ek,t& = 67.4 keV, 
plus the energy of the gamma, 6.2 keV, adding up 
to 73.7 keV. 

Lines at these energies can be seen in all spec- 
tra. For the crystal of which the spectrum is shown 
in figure 4, a measurement with a similar exposure 
to background radiation, about three months ear- 
lier, is available [8]. This allows us to extract the 
179 Ta and 181 W activities from two measurements. 
A comparison to the expected decay times is shown 
in figure 7, confirming the origin of these lines. 
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Figure 7: Two measurements of the 181 W and 179 Ta activi- 
ties in one crystal during two runs with similar exposure to 
background radiation. The expectations from decays with 
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and seen to be consistent with the measurement. Error bars 
are uncertainties of the fit to the spectral lines (each contain- 
ing tens of counts) and the running time of the respective 
measurements (detector Daisy, run 27 and run 28). 



6.5. Activated Calcium 

Capture of thermal neutrons can activate calcium 
present in the CaWC>4 crystals. 41 Ca has a half- 
life of 10 5 years and decays via electron-capture to 



41 K. The K-shell binding energy of potassium is 
only 3.61 keV, but our detectors are capable of de- 
tecting such low energies. Figure 8 shows again 
the spectrum of figure 5 but for lower energies and 
with a finer binning. We observe a prominent line 
at (3.71 ± 0.02) keV (mean and error from the fit of 
a Gaussian). This is only 3% above the expected 
energy from 41 Ca, which we therefore attribute to 
this line. The width of this line is only (130±19) eV, 
consistent with the calibration. For illustration, fig- 
ure 9 shows a typical pulse at this low energy, still 
clearly visible above baseline noise. 




Figure 8: Low-energy part of the spectrum shown in figure 5. 
In addition to the lines from 41 Ca and copper fluorescence 
there is a line at ~ 11.5 keV which so far could not be iden- 
tified in a consistent way (detector Verena, run 30). 
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Figure 9: A typical pulse with an energy of 3.7 keV (grey), 
and its 50-sample (2 ms) moving average (black). 

From the observed activity of (26 ± 4) /iBq we 
can readily calculate that there are (1.4±0.2) x 10 8 

Ca atoms present in the crystal, only a frac- 
tion of (2.2 ± 0.3) x 10~ 16 of all calcium atoms. 
This is an order of magnitude more accurate than 
other measurements of this fraction which are typ- 
ically carried out by accelerator mass spectrometry 
(AMS) [9, 10, 11]. Our present sensitivity is in the 
10~ 17 range even with CaW04 crystals that are not 
optimized for this purpose. This demonstrates the 
extreme sensitivity of our detectors for rare decays 
and opens new possibilities for radioactive dating 



using 41 Ca [12]. In fact, having an in situ cali- 
bration line at these low energies is an advantage 
for the direct Dark Matter search, since it indepen- 
dently confirms our energy calibration procedure in 
this otherwise inaccessible region. 

6.6. Copper Fluorescence 

Another line, seen at (8.17±0.04) keV in figure 8, 
can be attributed to copper fluorescence. This orig- 
inates from gammas or electrons that eject a pho- 
toelectron from the copper of the cryostat that sur- 
rounds the detectors, resulting in a copper X-ray 
absorbed in our crystals. Our energy calibration 
puts the line only 2% higher than the expected 
8.04 keV from the literature. This is consistent with 
the slight energy overestimate seen when identifying 
the 41 Ca spectral line, again confirming its origin. 

6.7. Lutetium 

All single crystal spectra above 100 keV are 
rather featureless. Yet, two crystals were operated 
simultaneously and in close vicinity during the com- 
missioning of the setup. This allows us to search for 
coincident events, of which a spectrum is shown in 
figure 10. 
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Figure 10: Spectrum of coincident events. To increase the 
counting statistics, the spectra of the two individual detec- 
tors were added to produce this figure. Without the coinci- 
dence requirement, the added spectrum would have a count 
rate of O(500) counts per 2keV bin (detectors Verena and 
ZORA, run 30). 

At 59keV, tungsten K Q escape events can be 
seen. In addition, two lines at 202 keV and 307 keV 
become very prominent, and we attribute those to a 
contamination with lutetium intrinsic to the crys- 
tals. 176 Lu beta-decays to 176 Hf with a total de- 
cay energy of 1191.7 keV. Each decay produces a 
gamma cascade with gamma energies of 306.8 keV, 
201.8 keV and 88.3 keV. In particular the higher 
energy gammas may escape from one crystal and 



hit the other. Hence, for a coincident 306.8 keV 
event observed in one detector, the other crystal 
will mostly absorb the other two gammas, plus the 
energy of the beta electron. This explains the struc- 
ture seen in the scatter plot of figure 11, show- 
ing the energies of coincident events in each de- 
tector. We observe a beta spectrum starting at 
201.8 keV + 88.3 keV = 290.1 keV for most events 
at 306.8 keV in the other detector (lines marked (a) 
and (b) in figure 11). Figure 12 shows the energy 
spectrum of these events in one detector, tagged by 
an energy of 306.8 keV in the other one, in agree- 
ment with the expectation from the lutetium beta 
spectrum. 



1(a) 




JtL"..'^!^.^....^^.-...'...!...:....!...'..-.^...! Uw 

100 200 300 400 500 600 700 800 
Energy / keV (First Detector) 



Figure 11: Scatter plot of coincident events in two detectors 
operated during commissioning. The vertical line (a) going 
up from 290 keV is due to lutetium decaying in the second 
detector, in coincidence with the 307 keV gamma escaped to 
the first detector. The horizontal line (b) is due to the same 
process with detectors reversed, at 307 keV in the second 
detector and starting at 290 keV in the first one. At (c) we 
observe a hint of a horizontal line starting at 395 keV due to 
lutetium decaying in the first detector with the energy of an 
escaped 202 keV gamma in the second detector. 



7. Conclusions 

The excellent energy resolution of less than 1 keV 
over the full energy range used in the Dark Matter 
search [1, 2] allows us to identify a variety of dif- 
ferent background sources in the CRESST experi- 
ment. Since the CRESST experiment can discrimi- 
nate nuclear recoil signal from electron recoil back- 
ground based on the light output, the presence of 
these features does not compromise the search for 
Dark Matter, yet allows an audit of the performance 
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Figure 12: Histogram of all coincident events in one detector 
that appear at (307 ± 3) keV in the respective other detec- 
tor: The 176 Lu beta spectrum starting at 290 keV is visible 
despite the low statistics. 
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of the experiment. The energy resolution (la) is 
130 eV at 3.6 keV, and the observed sharpness of 
the lines proves the high stability of the detectors 
during the background measurements running for 
several months. The observed line from activated 
calcium provides an independent in situ energy cal- 
ibration even below the energy range used for the 
Dark Matter search, and in addition proves the high 
sensitivity of the experiment to rare decays. 
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